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ABSTRACT 

rl [ Phases of elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix can 

be obtained using decays of B mesons to tt+tt", n^K^, and 7^~^K^ or ir'K . 
^ ; For B^ or 5^ — ^ 7r"'"7r~, one identifies the flavor of the neutral B meson 

;h I at time of production and studies the time-dependence of the decay rate. 

^ The other processes are self-tagging and only their rates need be measured. 

By assuming flavor SU(3) symmetry and first-order SU(3) breaking, one 
can separately determine the phases 7 = Arg and a = vr — /5 — 7, 
where (3 = Arg Vf^. Special cases include the vanishing of strong interaction 
phase differences between amplitudes, the possibility of recovering partial 
information when tt+tt^ and n^K^ decays cannot be distinguished from one 
another, and the use of a correlation between 7 and a in the region of allowed 
parameters. 
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I. INTRODUCTION 



The decays of B mesons [|I|] offer the prospect of confirming or refuting the current 
explanation of CP violation in the neutral kaon system p[, based on phases in the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix For example, unequal time-integrated 
rates for the tt+tt" decays of states which are initially and and B would signify CP 
violation, providing approximate information on the angle a of the triangle describing 
the unitarity of the CKM matrix. 

The presence of gluonic 0] and electroweak P penguin contributions in additional 
to the dominant ("tree") processes requires that one separate out several terms. An 
isospin analysis involving the study of the time dependence of the 'ir~^n~ mode and 
rates for the tt^h^ and n^n^ modes of B mesons, permits one to isolate the amplitudes 
contributing to final states with isospin and 2 and thereby to determine a rather well 
P, 0. However, for certain types of detectors, the observation of neutral pions may pose 
a challenge, and model calculations ^ predict a branching ratio for B^ of order 

10~^ or less. 

A few alternative ways to sort out the effects of several amplitudes in B^ tt+tt^ 
were suggested recently. DeJongh and Sphicas studied the dependence of the asym- 
metry in B^{t) TT~^7T~ on the magnitude and relative phase of the contributing terms. 
Using flavor SU(3) symmetry, Silva and Wolfenstein |]T0[ estimated the penguin con- 
tribution by comparing the tree-dominated decay rate of 5° tt^tt^ with that of 
B^ —* TT^K^ which has a large penguin term. Buras and Fleischer ]TT[ proposed re- 
lating the penguin term in B^ tt+tt^ via SU(3) to the time-dependent asymmetry 
of B^{t) K^K , where the penguin amplitude dominates. Kramer, Palmer, and Wu 



T^l note that the ratio of penguin to tree matrix elements is less model-dependent than 



either quantity alone, and thereby obtain a relation for a. Aleksan et al. [0] use model- 



dependent assumptions to learn the magnitude of the penguin effect on the measurement 
of a by relating the three AS* = decay modes tttt, Trp, and pp to the corresponding 
AS = 1 modes tiK, nK* and pK, and pK*. 

In this paper we examine in more detail a method proposed in Ref. il^ to determine 



phases of CKM matrix elements by detecting only kaons and charged pions in B meson 
decays. In the decays B^ — > tt^tt^ and B to vr^vr", one identifies the flavor of the 
neutral B meson at time of production and studies the time-dependence of the decay 
rate. One obtains the necessary information on additional amplitudes from the rates 
r(5° ^ n'K+), T(B° n+K-), and T{B+ n+K^) or 1(5" ^ n+lf) using flavor 



SU(3) symmetry [IIO], O, [TO, O, |I8| and first-order SU(3) breaking [jl9|. In the most 



general case we obtain information not only on a, but also on 7 = Kigiy*^ and on 
strong phase shift differences. Other ways to measure 7, based on charged B decays. 



were proposed in Ref. |20 



In Section II we describe the processes to be measured and the amplitudes on which 
they depend. We then study the precision to which various quantities can be determined. 
It is possible that the strong-interaction phase shift difference 5 between amplitudes is 
below detectable levels, in which case simplified analyses become necessary. Several of 
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these cases are discussed in Sec. III. The most general error analysis (for 5 7^ 0) is 
performed in Sec. IV, while Sec. V concludes. An Appendix is devoted to an aspect of 
Monte Carlo programs. 

II. PROCESSES AND AMPLITUDES 

A. Expressions for amplitudes and quantities quadratic in them 



We review the method proposed in Ref. which may be consulted for details. 
Our method employs flavor SU(3) symmetry [l^, |T^, and neglects "annihilation" 



amplitudes in which the spectator quark (the light quark accompanying the b in the 
initial meson) enters into the decay Hamiltonian |]18|. These amplitudes in B decays 
are expected to be suppressed by fs/^B, where fs — 180 MeV. We include first-order 
SU(3) breaking terms [1^, expected to be at most tens of percent, but neglect corrections 



expected to arise at a level of a few percent. 

In the SU(3) limit and neglecting annihilation terms, all B decay amplitudes into 
TTTT , ttK and KK states can be decomposed in terms of three independent amplitudes 
[]^, |18[: a "tree" term t(t'), a "color-suppressed" term c(c') and a "penguin" term p{p') 



These amplitudes contain both the leading-order and electroweak penguin p contribu- 
tions: 

t = T+{Cu- Cd)PEW ' C = C + (C„ - Cd)PEW , p = P + CdPsw ■ (1) 



Here the capital letters denote the leading-order contributions defined in Ref. |TB[, and 
Pew and Pew color-favored and color-suppressed electroweak penguin amplitudes 
defined in Ref. [0]. The values Cu = 2/3 and Cd = —1/3 are those which would follow 
if the electroweak penguin coupled to quarks in a manner proportional to their charges. 
(Small corrections, which we shall ignore and which do not affect our analysis, arise from 
axial-vector Z couplings and from WW box diagrams.) The AS* = amplitudes are 
denoted by unprimed quantities and the AS* = 1 processes by primed quantities. 

The amplitudes of the two processes 5° tt+tt^ and B^ tt^K~^ are expressed as 

A^^ = A{B''^n^n-) = -t-p=-T-P-'^P^^ , 

A^K ^ AiB' ^ n-K^) = -t' -p' = -r -P'- ^P'£^ , (2) 
while that for B^ —* n^K^ will be approximated by 

A^ ^ A{B+ TT+K') =p' = P'- ^-P'e%. ^P' + ^P'eSv , (3) 

neglecting a color-suppressed electroweak penguin effect of order |P^/P'| = 0((l/5)^) 
0. With this approximation, A^ contains the same combination of electroweak and 
gluonic penguins as in the expression for A^^k. 

The terms on the right-hand-sides of (0) and (|^) carry well-defined weak phases. 
The weak phase of T is ATg{VudV*i,) = 7, and that of P -|- \Pew approximately 
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ATg{VtdVtl) = — /?, where we neglect corrections due to quarks other than the top quark. 
The effects of the u and c quarks become appreciable when Vtd obtains its currently 
allowed smallest values. This corresponds to a small deviation of the CP asymmetry in 
B^{t) — » vr+vr" from sin(2a) sin(Amt) (where Am is the neutral B mass- difference). For 
large values of Vtd, where the deviation due to the penguin amplitude becomes significant 



22| , the u and c contributions become very small. T' also carries the phase 7, while the 



weak phase of P' + \P'ew is kig{VtsVg) = vr. 

In what follows we shall denote T = \T\, V = \P + ^PEwI^ ^' = \T'\, V = \P' + 
^Pew\- The ratio of AS = 1 to AS = tree and penguin amplitudes are given by 
the corresponding ratios of CKM factors, T'/T = |Vus/\4cz| = = 0.23, V'/V = 
\Vts/Vtd\ = Tf. To introduce first-order SU(3) breaking corrections, we note that in the 
\T'\ amplitude the W turns into an s quark instead of a in T. Assuming factorization 
for T, which is supported by experiments P3| , and justified for B ^ titt and tcK by 
the high momentum with which the two color-singlet mesons separate from one another, 
SU(3) breaking is given by the K/n ratio of decay constants 



T 

r 



1 


fK 


\Vud 





(4) 



Apart from small electroweak penguin terms, all amplitudes we consider are free 
of color-suppressed contributions, for which factorization might be more questionable. 
The situation would be very different were we to consider the amplitude for B^ tt^tt^, 
where the color-suppressed contribution could be dominant. 

In the penguin amplitudes (including electroweak penguin) of both B^ n^K^ 
and B^ ir'^K^ the b quark turns into an s quark instead of a. d in B^ —>■ tt+tt". 
Here we will denote the magnitude of the AS = 1 penguin amplitude by rtV, to allow 
for SU(3) breaking. Since factorization is questionable for penguin amplitudes, one 
generally expects V 7^ {fx/ f-K)^. We will assume that the phase 5p is unaffected by 
SU(3) breaking. Since this phase is likely to be small this assumption is not expected 
to introduce a significant uncertaintly in the determination of the weak phases. 

Assigning SU(3)-symmetric strong phases 5^, 5p to terms with specific weak phases, 
and taking account of SU(3) breaking, Eqs. (H) and (|]) may be transcribed as 

A^K = r„re^^-e*^ - nVe''^ , 

A+ = nVe'^p . (5) 

It will be shown that the numerous a priori unknown parameters in (|^), including the 
two weak phases a = vr — /5 — 7 and 7, can be determined from the rate measurements 
of the above three processes and their charge-conjugates. 

The amplitudes for the corresponding charge-conjugate decay processes are simply 
obtained by changing the signs of the weak phases 7 and (3. We denote the charge- 
conjugate amplitudes corresponding to (||) by A^^, A^^k, respectively. A state 
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initially tagged as a. oi B will be called B^{t) or B it). The time-dependent decay 
rates of these states to tt+tt" are given by 

+Im(e2^^^,,]4;jsin(Amt)] , 
n^it) ^ TT+TT-) = e-"[|A..|^sin^(^t) + |Z..|^cos^(^t) 

- Im(e2^^A^,A;^ sin(Amt)] . (6) 

Measurement of these quantities determines I^ttttP and \m.{e^^^ At^t^jC^^) . It is 

convenient to define sums and differences of the first two quantities, and we find 

A = + IZ^^I^) =r^ + 7'^-2rPcos5cosQ; , 

B = - \A^J^) = -2TVsm6sma , 

C = Im {e'^'f^A^^Al^) = -T^ sin 2a + 2TV cos 5 sin a , (7) 

where we use (3 + ^ — tt — a and where we define S = 5t — Sp. 

The rates of the self-tagging modes tt'K'^, ti^K~ and tt^K^ or 7r~K determine 
lA^rKp, IAttE'P and l^+l^, respectively. Again, we take sums and differences of the first 
two, and find 

D = 1-{\A^k\^ + \A^k\^) = {f^rf + V'^ - 2fuTV'cosScosj , 
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E = 1{\A^k\^ - \A^k\^) = 2f„r:P'sin5sin7 , 



= = :P'' . (8) 



2 



The rates for S+ — > tt'^K^ and B tt K are expected to be equal, since only 
penguin amplitudes are expected to contribute to these processes. Here we have defined 
P' = nV. 

Measurement of the six quantitities A — F suffices to determine all six parameters 
a, 7, T, V, V, 6 up to discrete ambiguities. The CKM parameter rt = \Vts/Vtd\, which 
is still largely unknown, is obtained from the unitarity triangle in terms of a and 7: 

sino; 

Tun^-. — . 9 

sm7 

We note that 

\a,k\' - \a^k\' = -(^)(|)(i^..r - i^^^n , (10) 

Jtt > 
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which determines the magnitude of SU(3) breaking in the penguin amphtude, V /V. 
The relation (^) between the particle-antiparticle rate differences in i? ^ t^K and in 
B ^ nn was recently derived p6[ in the SU(3) limit, fx/fn — * 1, V/V — 1. The 



authors assumed for SU(3) breaking a value V/V = Jk/ fn (based on factorization of 
penguin amplitudes) which in our approach is a free parameter to be determined by 
experiment. We expect it to differ from one by up to 30%. 



Both sides of Eq. (|T0|) are proportional to sin 6, and thus would vanish in the absence 
of a strong phase difference. In that case, one would have to assume a relation between 
V and V or some other constraint in order to obtain a solution. If, on the other hand, 

5 7^ 0, leading to a rate asymmetry between the self-tagging decays it~K^ and 
B — >■ n~^K~, the present method permits one to interpret that rate asymmetry in a 
manner independent of S. 

B. Likely ranges of observables 

The amplitude for B^ tt+tt^ is expected to be dominated by the T contribution, 
while those for B^ n^K^ and B^ tt^K^ are expected to be dominated by V. We 
shall choose units in which a branching ratio of 10~^ corresponds to a value of 1 for the 
rates A, D, and F. The normalizations of the other quantities are set accordingly. We 
shall also define the quantity 

A combined sample of the decays B^ — ^ tt+tt" and 5" tt'K^ has been observed 
with a joint branching ratio of (1.8l|]:^t[}J ± 0.2) x 10"^ [0, so that ^ = 1.8 ± 0.65. 
Equal mixtures of the two modes are most likely, corresponding to individual branching 
ratios of about 10~^ for B^ tt+tt^ and B^ —* tt^K^. A similar branching ratio is 
expected for B^ n'^K^ if the V amplitude dominates B ttK decays, as seems 
likely. Thus, values of order 1 for A, D, and F are expected. We shall consider a range 
of values for these quantities, subject only to constraints on the lower and upper limits 
for S. As we shall see in Sec. Ill D, when 6 can be neglected, this works out to a rule of 
thumb that + P'^ ^ 2. 

The detection of B^ n^K^ or B^ n^K^ in practice will utilize the channels 
B^ —>■ Tc^Ks, Ks TT+TT^, with a corresponding loss in efficiency of a factor of 3. We 
shall take this factor into account in estimating statistical errors on F. 

The remaining quantities B, C, and E are harder to anticipate. The Schwarz in- 
equality limits the value of |C| to be less than or equal to A. In practice we find values 
of \C\ larger than 2 to be very unlikely. Thus, we shall consider values subject to this 
restriction. Both B and E will vanish if 5 = 0. While a recent calculation |^8|] based 
on perturbative QCD p5| suggests that 5t 0, 6p ^ 9.5°, 6 ^ —9.5°, the possibility of 
non-perturbative effects (such as strong final-state interactions differing in channels of 
different isospin) cannot be excluded. Thus, we shall consider the representative values 

6 = 0, 5.7°, 36.9°, 84.3°, 95.7°, 143.1°, 174.3°. We take only non-negative values since 
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the error estimates are not affected by sign clianges in 6. The nonzero vafues wiff be 
discussed in Sec. IV. 



C. Constraints on the angles a and 7 



Recent analyses of constraints on the CKM parameters include those in Refs. ||29 
and [^]. We shall visualize the allowed ranges of a and 7 in order to choose illustrative 
sets of parameters. 



We begin with the Wolfenstein parameterization of CKM elements: 

V,t, = AX' , V^t, = AX%p~t7]) , Vt, = AX^{l-p-ir]) , (12) 
as well as others not quoted exphcitly, where A = 0.22. We shall assume fSO] Vcb = 



0.038 ± 0.003. (A slightly higher value is quoted in Ref. ||3^.) The measurement ||3^ 
\Vub/Vcb\ = 0.08 ± 0.02 based on charmless B decays implies (p^ + ry2)i/2 ^ g gg ^ o.09. 
The measurement oi B^ — B mixing implies \ytd\ = 0.009 ± 0.003, which we shall 
interpret as implying \1 — p — iri\ = 1.0±0.3. The requirement that the imaginary part 
of the — K mixing amplitude due to CKM phases be responsible for the observed 
CP violation in the kaon system implies a hyperbola 77(1 — p + 0.35) = 0.48 ± 0.20, 
where the 1 — p term in parentheses refers to the contribution of the top quark loop, 
while 0.35 refers to the charmed quark's contribution. 

The allowed region in (p, ■?]) is shown in Fig. 1(a); the corresponding range of (0,7) 
is depicted in Fig. 1(b) by a rather narrow band. The strong anticorrelation between 
a and 7 is a function of the limited range oi (3 = tt — a — 7, which is restricted to 
6.6° < P < 27° for the present set of parameters [0. The band would not be quite so 
narrow using the parameters of one other analysis ||29|| . 



Three representative points, noted on the figure, are described in Table I. These 
correspond to extreme and central values of a and 7. A number of illustrative examples 
will be presented for these points. 

Also shown on Fig. 1(b) is a linear least-squares fit to the points pi, p2, and p^: 
7 = 175° — 1.16a. Almost equally good is the approximate relation 

7 = 180° - 1.2a , (13) 

which we shall use in Sec. Ill A to simplify relations between vrvr and ttK rates. 

D. Limitation associated with size of 6 

The relation (|10D can be written as 

^ = -fjit . (14, 

We wish to evaluate the ratio V jV to better than 30% (the anticipated magnitude 
of SU(3) breaking). In this subsection we estimate the number of vr+vr" events (and. 
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Table I: Representative points in the (p, r]) plane and corresponding angles of the uni- 
tarity triangle. 



Point 


P 


V 


a 


P 


7 


n 








(deg.) 


(deg.) 


(deg.) 




Pi 


-0.30 


0.15 


20.0 


6.6 


153.3 


3.36 


P2 





0.35 


70.7 


19.3 


90.0 


4.16 


P3 


0.36 


0.27 


120.3 


22.9 


36.9 


6.35 



correspondingly, the size of the other data samples) needed in order to measure B and 
E [see the expressions (|^) and (^] to the required accuracy. 

Define the number of events averaged between particle and antiparticle decays: 



N{B'^ 71 + 71-) + iV(5° ^ 71 + 71- 



N(BO + N(B ^ ^^^j 

With equal branching ratios for tt+tt^ and 7T+K+ ^ the present data sample would consist 



of about 10 events each for A^^^^ and N^^k [271. The errors on A and B both scale as 



A'^/^, while those on D and E scale as N'Ik- Then in the samples of events used to 
measure A, B, D, and E, we expect 

6NA^6NB^Nli^ , 6Nd^6Ne-n'Jk ■ (16) 

We take as illustrative parameters T = V = V = 1, neglecting SU(3) breaking in 
the ratio V /V. Recalling the expressions for B and E, we expect the numbers of events 
in the samples corresponding to these quantities to be 

iV^B = — 2A^7r7r'"Msin5sin7 , Ne = ^Nt^k^u sin S sin j , (17) 

where we have used (P). Consequently, the fractional errors on Nb and Ne are 

\SNb\ 1 



I I 2Nni^ I sin 5 sin 7 1 
\6Ne\ 1 



I I - 1/2 • (19) 

\Ne\ 2A^7r4 r„| sin5sin7| 

Thus the fractional error on the quotient Nb/Ne is the sum in quadrature of these two 
errors: 

= I f 1 + II] (20) 

\{Nb/Ne)\ 2Ari/V„|sin5sin7| V fl) ' 
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Demanding that this error be less than 30% as noted above and substituting the values 
of the constants, we find 

91 

N^. > . 2, ■ 2 ■ (21) 
sm^ sm^ 7 

This gives an idea of the data samples required to improve upon the assumption of no 
more than 30% SU(3) breaking in penguin amplitudes. More detailed estimates are 
postponed until Sec. IV. Meanwhile, we examine the special case in which the strong 
phase shift difference 5 vanishes. 

III. VANISHING STRONG PHASE SHIFT DIFFERENCE 

Both the TTTT parameter B and the ttK rate asymmetry parameter E vanish when 
the strong phase shift difference 6 is zero. In that case, however, one can no longer use 
the relation (|r^) to determine the ratio V/V. One has 4 observables {A, C, D, and 
F) to determine 5 parameters (e.g., T, V, a, V, and 7). One must make additional 
assumptions to obtain solutions. In this section we explore several such possibilities. 

A. Simplified observables with 6 = 0. 

When 6 = 0, the equations in and (||) for A, C, and D become 

A = T^ + -2TVcosa , 

C = -T^sm{2a) + 2rVsma , 

D = {fuTf + P'^ - 2fuTV' cos 7 . (22) 
A simple relation follows from eliminating V between the first two of these equations: 

{C / Af = Az{l - z) , where z = r^sin^a/A . (23) 

The Schwarz inequality bound |C/A| < 1 mentioned earlier is manifest here. 

B. Linear relation between 7 and a 

A considerable simplification useful for anticipating the precision in determining a 
and 7 is obtained by noting that a and 7 are rather tightly correlated with one another 
[see Fig. 1(b)] as a result of the restricted range of (3. As mentioned, the dependence 
can be approximated by a straight line. If we are concerned mainly with learning the 
sign of p and are not so concerned about the exact magnitude of (p^ + •rfy/'^, we can 
substitute for 7 in the expression ( ^2]) for D, having already substituted F = "P'^, and 
thus each measurement of D implies a relation between T and a. 

An even greater simplification can be obtained if we neglect the term quadratic in T 
in D, and eliminate T, P, and V from the remaining equations involving A, C, D, and 



F. With the approximate formula (|T3D one has 



^"^(^•2^) = 2 6^ 1 (24) 

sm« Ja±7^^ 
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The sign ambiguity stems from the fact that the equation for T sin a has two solutions. 
We can anticipate that the solution with ^ A is the most likely, as long as V is 
relatively small compared to T, as generally anticipated. Since 

(rsina)2 = ^— , (25) 

and since \C\ tends to be small when sin a is near 1 (as for the point P2); "we anticipate 
that in that case we should choose the positive sign in the square root, and the argument 
of the overall square root in the denominator of the last fraction in Eq. (|2^) is about 
2A. On the other hand, when |C| is fairly large (e.g., for points pi and pa), the sign does 
not matter much, and the argument of the overall square root is about A. 

For 100 7r"''7r~ and 100 -r^K^ events, the errors in A and D are about 10%. Assuming 
that the \P'\ contribution is dominant in the error on F will then be about 17% 
(because of the branching ratio of neutral kaons to tt+tt^). One then finds an error of 
about 0.55 in the right-hand side of Eq. (p^ when a is near the middle of its range, 
and about 0.78 when a is near its lower or upper bounds. In Fig. 2(a) we plot the 
left-hand side of this equation, along with plotted points for a = 20° (pi), 71° (^2), 
and 120° (ps), with the errors in cos(1.2Q;)/sinQ; of ±(0.78, 0.55, 0.78), respectively. 
The allowed region in (a, 7) is shown in Fig. 2(b) along with the line corresponding to 
7 = 180° — 1.2a. The arrows designate values of a corresponding to pi, p2, and p^. 

When a is close to the center of its range, a sample of B decays corresponding to 
100 events in each of the 7i~^7i~ and n^K^ channels allows one to narrow the allowed 
region of a by roughly a factor of 2. For the lowest or highest allowed values of a one 
does somewhat better. For more precise estimates, one would retain the |Tp term in D 
when a ~ 90°, and would be more precise about the error on C. 

Let us for the moment neglect the small correction term in the expression (|^ for 
C. Then since sin 2a can take on the same value for two values of a equally above 
and below 7r/4, there is a discrete ambiguity associated with negative values of C and 
values of a < 7r/2. This ambiguity is likely to persist when we include the correction 
term. However, the addition of nK decay information appears capable of resolving 
this ambiguity, since it provides additional information on the angle 7 which is highly 
correlated with a. For values of a > 7r/2 and positive C only one solution (that for 
a < 37r/4) appears to be in the physical region, so we do not get the same sort of 
discrete ambiguity. 

C. SU(3) assumption for penguin amplitudes 

When (5 = 0, as mentioned, we are missing information on V'/V. In the previous 
subsection we supplied this information by assuming a functional relation between a and 
7. In the present subsection, we no longer assume such a relation, but simply assume 
this ratio to equal unity. (It was assumed to equal fx/ fn ~ 1-2 in Ref. |]10| and [p6[] .) 



Under this assumption, we may drop the tilde symbols on V in Eqs. (0). We may then 
plot contours of observables in the plane of T = |T| vs. V = \P'\ for various regions in 
the allowed parameter space. 
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Contours of fixed A are mainly sensitive to T, wliile those of fixed D depend mainly 
on V'. The slopes of the contours reflect the presence of constructive or destructive 
interference between T and V, depending on the signs of cos a and cos 7. Because of 
the strong anticorrelation between a and 7 shown in Fig. 1, the contours of A and D are 
nearly perpendicular to one another for each of the three illustrated cases. This means 
that for each pair [a, 7), a measurement of A and D selects a point in the (T, V) plane 
with comparable errors on T and V if the values and errors of A and D are comparable 
to one another. 

A measurement of F = \T"\^ (the n^K^ or t^^ K branching ratio, in units of 10~^) 
must be consistent with the determination just made. Thus, a one-dimensional allowed 
set of points (0,7) is chosen by the combined measurements of A, D, and F. The 
degree to which this choice is unique depends on being able to observe the effect of 
T — V' interference in the measurement of since in the absence of the contribution 
from T one would have D = F. That is, the average n^K"^ and n'^K^ (or tt^K ) rates 
would be the same in the absence of the T contribution to the tt^K^ mode. 

Once one has selected values of T and V for a one-dimensional set of points in the 
(a, 7) plane, the value of C can be used to distinguish among those points. Positive 
values of C tend to be associated with negative values of sin (2a) and hence with values 
of a greater than 90°. Such values correspond to {p,ri) values lying inside a circle of 
radius 1/2 with center atp=l/2, 7] = 0. These parameters correspond to roughly the 
right-hand one-third of the allowed regions in Figs. 1(a) and 1(b). The parameter spaces 
of Fig. 1 are much more sensitive to positive values of C than to negative values. We 
shall see such behavior again when we come to discuss a further simplification in the 
next subsection. 

D. Information without n/K separation 

The fact that the contours of A and D are nearly perpendicular to one another and 
have similar spacings in T and V', respectively, suggests that contours of 

S = A + D = {1+ fl)r^ + (1 + r-^)V''^ - 2TV\r-^ cos a + fnCOS7) (26) 

may not depend very much on which set of allowed (a, 7) one chooses. 

This expectation is borne out in Fig. 3. The observation that 5* is roughly indepen- 
dent of (a, 7) follows from the anticorrelation of cos a and cos 7. The sum of the average 

and 5^ branching ratios to tt^tt^ and tt^/^^ leads to an approximate constraint on 
7"2 -p/2 j^Q^g]^}y independent of CKM parameters within the allowed range. This is 
fortunate, since the CLEO Collaboration measures precisely this sum, with only 



weak distinction at present between pions and kaons. (Improved particle identification 
at CLEO is foreseen in the future.) One may expect a similar measurement in some 
hadron collider experiments, such as the CDF Detector at Fermilab, unless specific steps 
are taken for particle identification. 

A combined measurement of + V'^ and F = V'^ now can be used to determine 
each parameter. The determination of a and 7 is now more simple-minded, albeit less 
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precise, than in the previous section. Since our determination of T and V is now 
independent of a and 7, we can use these parameters in the equation ( P^ ) for C to plot 
contours of fixed C in the (a, 7) plane. (The variation with 7 springs from the fact that 
^u^^t = sin a/ sin 7, as mentioned previously.) An example of such a contour is plotted in 
Fig. 4 for the representative values T = V' = 1.1. One sees a fair amount of uncertainty 
for values of C near —0.7, where the contour intersects the allowed (a, 7) region in a 
wide range of points. This behavior is related to the discrete ambiguity noted at the end 
of Sec. Ill B. However, contours of positive C cut the allowed ger angle 

and lead to a more highly constrained solution. 

The measurement of C in the absence of particle identification is possible since one 
is following the time- dependence of a decay rate in which one compares the decays of 
states which were initially and i? to a combination of final states. The oscillations 
from which C is to be extracted are expected to stem only from the 7r"^7r~ final state. 

The imposition of particle identification returns one to the situation of the previous 
section. As mentioned, to make efficient use of the information thus provided, one must 
be able to see the difference between Z), where there is a small T contribution, and F, 
where there is none. 



E. Results of a Monte Carlo simulation 

We have explored numerically the case of 5 = 0, V = V using a Monte Carlo 
program which generates events with a statistical spread in the variables A, C, D, and 
F appropriate to data samples corresponding to a total of M decays of or B to 
TT+vr". Scaling other quantities to the expected vr+vr^ rates, and recalling the discussion 
of Sec. II D, we then assume 

5A = 5C = 5D = 5F/V3 = l/^M . (27) 

For M = 100, 1000, 10000 we then ask how well a and 7 can be determined. Numerically 
this is accomphshed by stepping a and 7 through a range of values, accepting any 
solution which is within 2cr of the generated value for each parameter A, C, D, F, and 
averaging all such solutions. To allow for the possibility of multiple solutions, a cluster 
algorithm (described in the Appendix) is applied. We also restrict 10° < a < 130°, 
20° < 7 < 170°, and I7 — (175° — 1.16a) | < 30° in accord with the allowed regions in 
Fig. 1. The results are shown in Fig. 5. 

One sees a noticeable improvement with increased statistics. A clear distinction 
between the cases (pi or P2) and ps is already possible with M = 100. Satisfactory results 
for all three points are obtained for M = 1000. The small ambiguity for pi appears to 
be related to the multiple intersections of the contours for negative C (Fig. 4) with the 
allowed region of parameter space. Errors are reduced further when M is increased to 
10000. 
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IV. ERROR ANALYSIS FOR GENERAL FINAL-STATE PHASES 



In Ref. |T^ we estimated the statistical accuracy of determining the weak phases a 



and 7 using the present method to be at a level of ten percent, given around a hundred 
events in each channel. The theoretical uncertainty of the method is at a similar level, 
involving the following corrections all of which are of order a few percent: A correction 
from an electroweak penguin amplitude in n^K'^, corrections due to u and c 

quarks in the —>■ ti^ti~ penguin amplitude, second-order SU(3) breaking in the 
magnitudes of weak amplitudes, first order SU(3) breaking in the (small) strong phase 
of the penguin amplitude, and O^fsI'mB) annihilation amplitudes. In this subsection 
we analyze in more detail the precision on a and 7 that can be attained with a given 
sample of events as a function of the parameters. We use a Monte Carlo program similar 
to that described in Sec. Ill E to generate events with a given Gaussian distribution in 
the parameters A — F appropriate to a total sample corresponding to M ti^ti~ decays. 
We choose T = 1 and V = V = 1 for the purpose of generating events. In addition to 
the errors assumed in Eq. (^^, we assume 



5B = 5E = 1/VM . (28) 

The results are shown in Fig. 6. 

The method clearly improves with increased statistics, in a manner roughly compat- 
ible with our estimate (|2l|) . Despite the presence of a large spread in values of a and 
7, one can already distinguish the case ps from the cases pi and p2 for M = 100. The 
distinction between pi and p2 appears to emerge by the time one reaches M = 1000. 
On the other hand, one sees the distinct appearance of clusters of points, corresponding 
to the presence of discrete ambiguities. These ambiguities are to be distinguished from 
the one mentioned at the end of Sec. Ill B, and appear to be related to uncertainties 
in the value of the final-state phase. Their nature is best ascertained by referring to 
the case M = 10000. One can detect up to three clusters of solutions. For example, in 
the case of input parameters in the vicinity of pi, a final-state phase 6 ~ 90° turns out 
to be ambiguous with two other cases, one with 6 < 90° and the other with 6 > 90°, 
as one sees by referring to the corresponding plots of a vs. 6 shown in Fig. 7. For 
input parameters near p2, the more serious ambiguities appear to occur for moderate 
or small values of 6. For input parameters near p^, uncertainties present for M = 100 
and M = 1000 appear to have largely disappeared by the time M reaches 10000. This 
behavior may be related to the uniqueness of the solution provided by large positive C 
for points near p^, in Fig. 4, but also points to the absence of ambiguities associated with 
the value of 6. 

The discrete ambiguities mentioned previously are quite noticeable in Fig. 7; in 
addition, for S near 180°, even the largest value of M does not lead to solutions in which 
that value is uniquely determined. 
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V. CONCLUSIONS 



To summarize, we have shown that measurements of the rates for B decays to modes 
involving charged pions and kaons in the final states can determine the shape of the 
unitarity triangle, even in the absence of theoretical or experimental information about 
final-state phases. The full set of measurements involves the detection of the time- 
dependent rates for and B — > tt+tt' and the rates for B° n-K+, B n+K', 
and B^ Kstv^- A rate asymmetry between S° — > n'K'^ and B 7i~^K~ is 
needed in order to perform a solution for all necessary parameters. In the absence 
of this asymmetry, one can obtain partial information by noting the tight correlation 
between the angles a and 7 of the unitarity triangle, or by assuming an SU(3) relation 
between strangeness-changing and strangeness-preserving penguin amplitudes. One can 
even dispense with particle identification, summing tt+tt" and n^K^ modes, if only 
crude constraints on parameters are desired. As a result of the strong anticorrelation 
between a and 7 in the physically allowed region of parameter space, the nK modes are 
particularly helpful in resolving a discrete ambiguity associated with the behavior of the 
function sin 2a which would be present if one studied tttt modes alone. 

In the simplest case examined, where the assumption of 5 = and the strong corre- 
lation between a and 7 in the allowed parameter space were utilized, we found that a 
sample of events corresponding to 100 tt+tt^ and 100 -k^K^ events, with a correspond- 
ingly reduced number of detected B^ Ksti^ decays, was sufficient to reduce the 
allowed region in parameter space by roughly a factor of two, depending on the values 
of the CKM angles. 

In the more general case in which 5^0 but no relation between a and 7 was 
assumed, we found that by assuming SU(3) symmetry for penguin amplitudes we could 
obtain unique solutions for a and 7, with some possibility of discrete ambiguity when 
a is small and 7 is large (corresponding to p < in the language of the Wolfenstein 
parametrization). Even when a distinction between charged pions and charged kaons is 
not possible (Sec. Ill D), partial information on the parameters can be obtained, since 
the time-dependent effects are expected to be confined to the tt+tt^ channel and thus a 
measurement of the parameter C (defined in Sec. II) is still possible. 

In the most general case of nonzero final-state phase differences 5 we find that the 
program described here requires approximately 10^/(sin^ 5 sin^ 7) decays of neutral S's 
to charged pions (and a similar number of ttK events) in order to free oneself from 
assumptions of SU(3) breaking at the 30% level in penguin amplitudes. A Monte Carlo 
program has shown that one begins to get useful information with 100 such decays (to 
be compared with about 10 in the present data sample). The full power of the method 
becomes apparent as the sample exceeds 1000 and approaches 10000. Even under such 
circumstances, a discrete ambiguity remains associated with the size of final-state phases. 
Arguments external to those presented here (such as the allowed regions in the (a, 7) 
parameter space, the expected magnitude of SU(3) breaking, and the expected size of 
final-state phases) may be necessary to resolve such ambiguities. 
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APPENDIX: DETAILS OF CLUSTER ALGORITHM 

Given a set of values for A, B, C, D, E and F, the values of a, 7 and 6 are not neces- 
sarily determined uniquely. Apart from the ambiguities associated with the numerical 
nature of the algorithm, there can also be discrete ambiguities. In that case the set 
of triplets (a, 7, 5) consistent with all the observed quantities will form clusters in the 
(a, 7, 6) space, for any given set of A, B, C, D, E, F. The number of clusters correspond 
to the number of discrete solutions and the spread within a cluster corresponds to the 
(numerical) error on that particular point. The average of all points in each cluster is 
taken to be the central value for that cluster and is plotted in Figs. 5, 6 and 7 as a single 
point. The number of points plotted for each data set is thus the number of discrete 
solutions for that data set. 

The ambiguities associated with the numerical nature are expected to be continuous; 
i.e. for any point i to belong to a cluster, there should be at least one point j in the 
cluster such that \ai — aj| < Aa, where Aa is the least count in a in the numerical 
algorithm. It is observed that this condition alone is sufficient to separate the clusters 
and additional conditions on 7j or 6i are not needed. Two points i and j belonging 
to different clusters will fail to satisfy this condition. Different clusters can thus be 
separated from each other. 
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FIGURE CAPTIONS 



FIG. 1. Allowed ranges of CKM parameters (bounded by solid lines). The points pi, p2, 
and p3 are described in Table I. (a) (p, r]) plane; (b) as function of angles a and 7 of the 
unitarity triangle. A linear least-squares fit to pi, p2, and ps yields 7 = 175° — 1.16q;, 
as shown by the dashed straight line in (b); its map into the (p, ?]) plane is the dashed 
curve in (a). 

FIG. 2. (a) The function cos(1.2a)/ sin a as a function of a, and errors on a expected for 
a sample corresponding to 100 or — > 7r"'"7r^ decays. The three points, from left to 
right, correspond to pi, p2, and ps of Table I, and 5 = is assumed, (b) Corresponding 
regions in a — 7 plane. 

FIG. 3. Contours in the |T| — \P'\ plane of fixed S = A + D (sum of average tt~^7t~ 
and TT^X^ branching ratio in units of 10~^), for 6 — 0. Dotted curves: S — 0.2. Other 
curves, outward from origin: S — 1, 2, 3, 4. 

FIG. 4. Contours of fixed C for |T| = \P'\ = 1.1 as functions of a and 7, for 6 = 0. Only 
those branches intersecting the allowed region are shown. Dotted curve: C = —0.7; 
dashed hyperbola: C = —0.5; sohd curves, from left to right: C = 0, 1; dashed eUipse: 
C = 1.5. 

FIG. 5. Scatter plots in the a — 7 plane of 100 events generated according to errors 
appropriate to samples of M = 100, 1000, 10000 events for the points pi, p2, Ps of 
Table I, for (5 = 0. Here the a and 7 axes are plotted in degrees. 

FIG. 6. Scatter plots in the a — 7 plane for nonzero input values of 6 (labels above 
columns) for the points pi, p2, pa of Table I (labels to left of rows). Here 0° < (a, 7) < 
180°. (a) M = 100; (b) M = 1000; (c) M = 10000. 

FIG. 7. Scatter plots in the a — 6 plane for nonzero input values of 6 (labels above 
columns) for the points pi, p2, pa of Table I (labels to left of rows). Here 0° < a < 180°; 
-180° <6< 180°. (a) M = 100; (b) M = 1000; (c) M = 10000. 
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